CERN-TH/99-165 

On Average Properties of 
Inhomogeneous Fluids in General Relativity 

g^; I: Dust Cosmologies 

o\ 

o: by 

Q 

^ '. Thomas Buchert 



(N 

> 

in 



Theory Divison, CERN, CH-1211 Geneve 23, Switzerland * 



\^ , e-mail: bucliert@tlieorie.pliysik.uni-niuenchen.de 

O 

a^ 
o> 

^ I Slimiliary ! For general relativistic spacetimes filled with irrotational 'dust' a 

PfJ' generalized form of Priedmann's equations for an 'effective' expansion factor ax>(t) 

of inhomogeneous cosmologies is derived. Contrary to the standard Friedmann 
K^ ! equations, which hold for homogeneous-isotropic cosmologies, the new equations 

H I include the 'backreaction effect' of inhomogeneities on the average expansion of 

the model. A universal relation between 'backreaction' and average scalar cur- 
vature is also given. For cosmologies whose averaged spatial scalar curvature is 
proportional to a^ , the expansion law governing a generic domain can be found. 
However, as the general equations show, 'backreaction' acts as to produce aver- 
age curvature in the course of structure formation, even when starting with space 
sections that are spatially flat on average. 



* 



on leave from: Theoretische Physik, Ludwig-Maximilians-Universitat, Theresienstr. 
37, D-80333 Miinchen, Germany 





1. The Averaging Problem 



A longstanding question in cosmology is, how to average a general inhomogeneous model 
(Ellis 1984). Also one would like to know under which assumptions, if any, the average vari- 
ables obey Friedmann's equations lying at the basis of any theory for structure formation 
in the Universe. An answer for cosmologies containing a pressure-free fluid ('dust') has 
been given recently in the framework of Newtonian cosmology (Buchert & Ehlers 1997; see 
the references therein for the formulation of the averaging problem and earlier attempts to 
solve it). Some results relevant to the present work may be briefly summarized as follows: 
Consider any simply-connected, compact spatial domain in Euclidean space V C lE"^ with 
volume \T>\ = Vt>- Then, the spatial average of Raychaudhuri's equation for the evolution 
of the expansion rate, under the assumption of conservation of the domain's mass, yields 
an equation for the scale factor ax? oc V^y' (which depends on content, shape and posi- 
tion of the spatial domain). This equation contains as source terms, besides the average 
mass density, averages over fluctuations of the shear, vorticity and expansion scalars due 
to the presence of inhomogeneities. These 'backreaction' terms vanish, if the average is 
performed over the whole universe having topologically closed space sections. For van- 
ishing 'backreaction' this equation is equivalent to the standard Friedmann equation for 
homogeneous-isotropic universes. 

In the present paper I provide the corresponding answer in the framework of general 
relativity. A geometrical relation, having no Newtonian analogue, delivers additional in- 
formation that allows to obtain a single equation relating the 'backreaction' and average 
curvature terms. This equation holds for any spatial domain and for a large class of inho- 
mogeneous cosmologies with curved space sections without perturbative assumptions. The 
solution of the 'backreaction problem' for scalar characteristics can be found in the case 
of space sections whose Ricci scalar averages out to zero, or, displays a dependence oc a^^ 
as in the standard model, which implies that the averaged scalar curvature decouples from 
the 'backreaction' term. However, in general, 'backreaction' due to the presence of inho- 
mogeneities impacts on the averaged Ricci scalar in the course of structure formation. A 
general solution to this problem would incorporate a scale-dependent description of the 
density, the expansion and other scalar variables of any structure formation model. 

In Section 2 we proceed by prescribing the basic equations and the averaging pro- 
cedure; then we give the general equations governing the domain-dependent scale factor 
ax>{t) in a Theorem. An equivalent set of these equations presented in Corollary 1 demon- 
strates the surprising result that, inspite of the non-commutativity of averaging and time 
evolution, the averaged quantities obey the same equations as the local ones. Inferred from 
the averaged equations Corollary 2 defines a set of average characteristics in analogy to the 
cosmological parameters of the standard model. A compact form of the previous results is 
presented in Corollary 3, which displays a universal relation between average scalar cur- 
vature and 'backreaction'. Finally, in Section 3, we discuss some immediate consequences 
of this result, and comment on related work. In Appendix A we present the Newtonian 
analogues, while Appendix B gives some illustrative examples of expansion laws. Appendix 
C is dedicated to an alternative derivation of the averaged equations. 
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2. Averaging Einstein's Equations for Scalars 



2.1. Choice of foliation and basic equations 

We shall restrict ourselves to the case of irrotational fluid motion with the simplest matter 
model 'dust' (i.e. vanishing pressure). In this case the flow is geodesic and space-like 
hypersurfaces can be constructed that are flow-orthogonal at every spacetime event in a 
3+1 representation of Einstein's equations. 

We start with Einstein's equations^ 

R^^u - -^dfjiuR = StiGqu^u^ - Kg^^ , (la) 

with the Ricci tensor R^,^, its trace R, the fluid's 4— velocity u^ {u^u^ = —1), the cosmo- 
logical constant A, and the rest mass density g obeying the conservation law 

(Qu^un ,, = . (16) 

We choose a flow-orthogonal coordinate system x^ = {t^X'^) (i.e., Gaussian or normal 
coordinates which are comoving with the fluid). Writing x^ = /^(X^,t) we have -u^ = 
f^ = (1, 0, 0, 0) and u^ = f^ = (—1, 0, 0, 0), where the dot denotes partial derivative with 
respect to proper time t. 

These coordinates are deflned such as to label geodesies in spacetime, i.e., u'^u'^.^, = 0. 
With the choice of vanishing 3— velocity the coordinates are in addition chosen to be 
comoving. Thus, in a 3 + 1-splitting of spacetime, the spatial set of Gaussian coordinates 
also labels fluid elements or trajectories in 3-space, X^ = 0; we are entitled to call X'^ 
Lagrangian coordinates, because they are identical to those in classical fluid dynamics. It 
should be emphasized, however, that the flnal result will be covariant with respect to the 
given foliation of spacetime and thus not dependent on this particular choice of coordinates. 

Let (t, X^) be the independent variables. As dependent variables we may choose the 
spatial 3— metric gij (the flrst fundamental form of the hypersurfaces of constant t) in the 
line element 

c/s2 _ _^^2 ^ g^^dX'dX^ , (Ic) 

the extrinsic curvature tensor Kij := —h^^h -Ua-^p (the second fundamental form of the 
hypersurfaces of constant t) with the projection tensor into the hypersurfaces orthogonal 
to Wq,, haf3 := (7a/3 + ^aUp, and the rest mass density g. Einstein's equations (la) together 
with the continuity equation (lb) (contracted with Ui,) then are equivalent to the following 



^ Greek indices run through 0...3, while latin indices run through 1...3; summation over 
repeated indices is understood. A semicolon will denote covariant derivative with respect 
to the 4-metric with signature (— , +, +, +); the units are such that c = 1. 



system of equations (see, e.g., Arnowitt et al. 1962, York 1979) consisting of the constraint 
equations^ 



1 '— "2 _ T^i J^j 



-in + K^- K'^K\] = SnGg + A , (2a) 

^^•||-^b=0 , (26) 

and the evolution equations for the density and the two fundamental forms: 

g = Kg , (2c) 

{g.jY =-2 g,uK^^ , {2d) 

{K\) ' = KK'^ + TV^ - {AtxGq + A)5^ . (2e) 

IZ •.= 7V' ^ and K := K\ denote the traces of the spatial Ricci tensor TZij and the extrinsic 
curvature tensor Kij, respectively. Expressing the latter in terms of kinematical quantities. 
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-Kij = Oij = a^j + -Ogij ; -K = 6' , (3) 



with the expansion tensor G^j , the trace-free symmetric shear tensor aij and the expansion 
rate 9, we may write Eqs. (2) in the form 

\n + \e'^ -a'^ = %tiGq + K , (4a) 

Q = -eQ , (4c) 

{g,,y = 2 g.kcr^^ + -%fc5^j , (4rf) 



(a^)•=-^a^-7^v + ^5^ 



a'^ -\e'^ + ^tiGq + K 



(4e) 



where we have introduced the rate of shear a^ := |(T* a-'j. To derive this last equation we 
have used Raychaudhuri's equation 

e + \e'^ + 2a'^ + ATiGQ-K = , (5a) 

which follows from the trace of Eq. (4e) combined with the constraint (4a) . Using again 
the constraint (4a), we may cast the trace-free part (4e) into the form 

Kr+^^^• = -fe-^5^7^) . (56) 



^In Eq. (2b) n denotes covariant derivative with respect to the 3— metric gij, while a 
single vertical slash denotes partial derivative with respect to the Lagrangian coordinates 
X*; note that in the present case the covariant spatial derivative of a scalar reduces to the 
partial derivative. The overdot denotes partial time-derivative (at constant X*) as before. 



This set of equations has recently been discussed in connection with perturbation theory 
by Kasai (1995), Matarrese (1996, and ref. therein) and by Matarrese & Terranova (1996), 
as well as in the papers by Russ et al. (1996, 1997). Here, we proceed without perturbation 
theory. 

Taking the trace of Eq. (2d), written in the form 
and defining 



J(t,X^):=^det((7.,) , (6a) 

we obtain with ^(7*'^ [gki]' = (In J)' the identity 

j = -KJ = ej . (66) 

Using it we can integrate the continuity equation for the rest mass density (2c) along the 
flow lines: 

p(t,X^) = (^(to,X^)J(to,X^))J-^ . (7) 

Both, Raychaudhuri's equation (5a) and the integral of the continuity equation (7) are 
identical to their Newtonian counterparts (Buchert & Ehlers 1997). Below we shall make 
explicit use of the constraint (4a) , which has no Newtonian analogue. This equation will 
provide a key element for the understanding of the 'backreaction' problem. 

For later discussion it is convenient to also introduce the abbreviations I and II for two of 
the scalar invariants of the expansion tensor, its trace, 

l-=@^^ = e , (8a) 

and the dispersion of its diagonal components. 



2.2. Averaging the traces of Einstein's equations 

Spatially averaging equations for scalar fields is a covariant operation given a foliation of 
spacetinie. Therefore, as we already pointed out in the paper on averaging the Newtonian 
equations (Buchert & Ehlers 1997), we may average, e.g., Raychaudhuri's equation (5a) in 
full formal analogy to the Newtonian case. 

Let us define the averaging operation. Spatial averaging of a scalar field \1/ as a 
function of Lagrangian coordinates and time on an arbitrary compact portion of the fluid 
V is straighforward"^ and is defined by the volume integral 

(vl>(t,XO)x) := 7^ / Jd^X vl>(t,X^) , (9a) 

with the volume element dV := Jd^X of the spatial hypersurfaces of constant time. The 
volume itself is given by^ 

Vv{t) := / Jd^X . (96) 

Jv 

We also introduce a dimensionless ('effective') scale factor via the volume (normalized by 
the volume of the initial domain Vt)^ ) , 

Thus, the averaged expansion rate may be written in terms of the scale factor: 

Vx> (iv 

The integral (7) states the conservation of the total rest mass Mx> within a portion of the 
fluid T) as it is transported along the flow lines. 



Mjy = / Jd'^X Q = const. ^ (g)T) = ^ . (9e) 



Employing this averaging procedure we may easily prove many statements found in 
(Buchert & Ehlers 1997) which also hold in general relativity. From these results we 
are going to use the 'Commutation rule' (here written for an arbitrary scalar fleld \1/): 



"^This averaging method functionally depends on content, shape and position of the 
spatial domain of averaging, which we consider as being given (see: Stoeger et al. 1999 
for alternative averagers). It should be stressed that we do not attempt to average the 
spacetime geometry as a whole; useful information for cosmology may be already obtained 
by averaging the scalar parts. 

^In comparison with the Newtonian deflnition Vx> = Jx>(t) ^^^ ~ Iv(t ) Jd^X, where 
Xi are Eulerian and Xi Lagrangian coordinates, our domain V corresponds to a Lagrangian 
domain, because it is transported along geodesies; however, in contrast to the Newtonian 
case, it is implicitly time-dependent due to the evolution of the 3-metric. 



Lemma (Commutation rule) 

{^)h - {^)v = {^e)v - {^)v{e)v . (9/) 

Averaging the Hamiltonian constraint (4a) and Raychaudhuri's equation (5a) with the 
help of the prescribed procedure, we end up with the following two equations, which we 
may summarize in the form of a theorem. 

Theorem (Equations for the effective scale factor) 

The spatially averaged equations for the scale factor a-D, respecting mass conservation, 
read: 



averaged Raychaudhuri equation: 



■D 



3— +47rG-^^-A = Qt, ; (10a) 



averaged Hamiltonian constraint: 



^T)\ „ ^ M-D 1,^\ . Qv 



3^ _87rG— ^ + -(7^)x,-A = -^ , (106) 

where the mass Mj), the averaged spatial Ricci scalar {TZ)ti and the 'backreaction' Q-p are 
domain-dependent spatial constants and, except the mass, time-dependent functions. In 
particular, the 'backreaction' source term is given by 

Qx):=2(II)x)-^(I)^ = ^((^-(^M')x,-2(a% . (10c) 



We also note the following surprising property of the averaged equations compared with 
their local forms: in spite of the non-commutativity of the averaging procedure and the 
dynamical evolution, which is expressed by the Commutation rule (9f), we find that the 
same equations hold for the averaged and the local quantities provided we express them 
in terms of the invariants (8) . This establishes the following 

Corollary 1 (Averaged equations) 

The spatial averages of the Hamiltonian constraint (4a) , the continuity equation (4c) and 
Raychaudhuri's equation (5a) read: 

^(7^)x) = ^nG{Q)T, + A - (11)^ , (11a) 

{Q)h = -{0)v{Q)v , (116) 

{e)^ = -{e)l + K-AnG{Q)T, + 2{ll)T, , (lie) 

i.e., the averages (^)d, (6')x), {T^)v and (II)i? obey the same equations as the local fields q, 
9, TZ and II. (The reason for this nontrivial property is the special type of nonlinearities 
featured by the gravitational system, e.g. the nonlinearity in 9 contained in Raychaudhuri's 
equation.) 
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Corollary 2 (Dimensionless characteristics of inhomoqeneous cosmologies) 

As in the standard homogeneous-isotropic cosmologies we may introduce a domain- 
dependent Hubble function Hxi := ^^, and dimensionless average characteristics as follows: 

Oa : = ^ , (126) 



'I? 



r> 



^fc ■■ = -TTTTf , (12c) 

which, in view of (10b), obey 

n^ + Oa + Ofe + Oq = 1 . (12e) 

All these dimensionless "cosmological parameters" actually depend on the spatial scale 
of averaging including the dimensionless cosmological constant, which depends on scale 
through H-i)- 

The equations (10a,b) form a system of two equations for the three unknown variables 
d-D, {T^)v and Q-d- Therefore, we cannot solve the 'backreaction' problem for scalars based 
on this system. We may eliminate the 'backreaction' term from (10a) and insert it into 
(10b). This results in an equation for the scale factor aj) and the average Ricci scalar of 
the domain. Alternatively we may proceed as follows: we calculate the time-derivative of 
Eq. (10b) and insert into the resulting equation (10a) and (10b). This yields a universal 
relation between the averaged Ricci scalar and the 'backreaction' term: 

Corollary 3 (Relation between average scalar curvature and 'backreaction') 

A necessary condition of integrability of Eq. (10a) to yield Eq. (10b) is provided by the 
relation: 

Qh + ^—Qv + {n)h + 2'^{n)T, = , (13a) 

av a-D 

or, equivalently, 

(a|,QI,)•+a|,(a|,(7^)x))• = . (136) 

In Appendix B we give some examples of expansion laws that can be derived from this 
relation. 



Notes: 

Eq. (10a) has been already confirmed as a valid equation in general relativity in the work 
on averaging the Newtonian equations (Buchert & Ehlers 1997); an equation analogous 
to Eq. (10b) has been derived by Carfora & Piotrkowska (1995) in connection with man- 
ifold deformations at one instant using the constraint equations of general relativity. In 
our derivation of Eq. (10b) we have inserted the 'backreaction' term Qti back into the 
constraint equation (4a) and have used (9d). 

Russ et al. (1997) have also used the equations (10a,b) (in a truncated form and 
using a reference background solution) for the purpose of perturbative calculations of 
the 'backreaction' effect. Note that, contrary to their derivation, we have not performed 
a conformal rescaling of the metric, nor have we used the splitting into a background 
reference solution and deviations thereof. We postpone further comments on their work 
to Section 3. 



3. Discussion and Perspectives 



3.1. Summary of results 

We have derived a generalized form of Friedmann's differential equations including 'backre- 
action terms' due to the presence of inhomogeneities. One of these equations was obtained 
on the basis of averaging Raychaudhuri's equation on spatial domains whose mass content 
is preserved in time. It is formally identical to the equation derived in the framework of 
Newtonian cosmology (Appendix A). The other equation arises by averaging the Hamilto- 
nian constraint (having no Newtonian analogue). It delivers an additional relation between 
the averaged spatial Ricci curvature scalar and the 'backreaction term'. We have shown 
that there exist special solutions which describe the evolution of the average curvature and 
'backreaction' terms (Appendix B). 

Let us list some immediate conclusions which may be drawn: 

1. The average expansion of inhomogeneous cosmologies does, in general, not follow the 
expansion law of the standard FRW cosmologies. There are, however, generalized expan- 
sion laws which govern the motion of arbitrary spatial domains, provided assumptions on 
the relation between the averaged Ricci scalar and the 'backreaction term' are made. Here, 
perturbation theory would be useful to establish such relations. 

2. The general expansion law shows that 'backreaction' due to the presence of inhomo- 
geneities implies the existence of a non-vanishing average Ricci scalar in general situations. 
This is true even if we consider domains which are on average Ricci fiat at some initial 
instant. 

Some comments about these conclusions are in order. 

adl. We may stipulate that the assumption of vanishing average Ricci scalar could be 
a sensible one, if we consider typical portions of the Universe (which itself may have on 
average flat space sections to a good approximation). As demonstrated in Appendix B a 



solution can be obtained in this case and the evolution of the average expansion is then 
exactly known. Consideration of the general case reveals, however, that this point of view 
is too naive: Looking at Eq. (10b) we must expect that the dimensionless contribution to 
the averaged scalar curvature {Ti)v/G{Q)v might be of the same order as the dimensionless 
contribution to the 'backreaction' Qv/G{Q)-r)- Approximating the average curvature by 
zero relies on a similar prejudice as saying that 'backreaction' may be neglected. Both 
approximations imply restrictions on general inhomogeneous models. 

ad2. Here, we may imagine the likely situation that an initially critical universe in the 
sense of the Einstein-de Sitter model may evolve into an under- or overcritical universe, 
respectively, in the course of structure formation. Hence, it is possible that a Ricci flat 
universe develops into a universe with on average negative/positive spatial curvature at the 
present epoch. From the point of view of the standard inflationary paradigm the former 
situation is favoured, when the theoretical expectation of an on average Ricci flat universe 
at the exit epoch is compared with measurements of the density parameter at the present 
epoch. 

For the general case a solution seems to lie beyond the scope of this article. Let us 
illustrate why the system of averaged equations (13a,b) cannot be closed on the level of 
scalars and also, how we may achieve closure by additional assumptions. 



3.2. Attempting closure of the averaged equations 

In order to obtain a more general result, we would like to find an independent evolution 
equation for the spatial Ricci scalar. Kasai (1995, appendix) has derived an evolution 
equation for the spatial Ricci tensor. His equation reads^: 

(7^^) • - 2K^7^^ = -K^iif - k] %, + K\ %, + < \ . (14a) 

This relation is purely geometrical and makes no use of the field equations. Taking the 
trace of Eq. (14a) we first obtain 

1^ .,^ _^ /.,£ 



-7^-K<^,7^^ = -(^K^^ll,-K|.j . (146) 

The r.-h.-s. of this equation vanishes according to the momentum constraints (2b). Using 
these constraints, the field equation (2e) to eliminate TZij in favour oiKij and the definition 
(8b) we obtain with Kij = —Qij'- 

-n = ee + e^-ii-2eii-e{A'KGQ + K) . (i4c) 

Inserting Raychaudhuri's equation (5a) we find that Eq. (14c) is just the time-derivative 
of the Hamiltonian constraint (4a) combined with (4c) . Hence, the trace of the evolution 



^Note that we have used the canonical definition of the extrinsic curvature in this 
paper: Kij = -&ij. 
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equation (14a) for the spatial Ricci tensor cannot be used to close the system of equations 
and to solve the 'backreaction' problem. 

We may try to use Eq. (5b), which is an equation for the trace-free parts, and contract 
this equation into a relation among scalar quantities. Indeed, if we contract Eq. (5b) with 
a^'j and eliminate the expression TZ^ .a-'^ from Eq. (14b) we obtain 

Using the Hamiltonian constraint (4a), Eq. (14d) can be written as 

7^ + ^7^ - 2M + 211 + 2^11 = , (14e) 

which, however, can also be obtained by inserting the Hamiltonian constraint into its 
time-derivative. 

These examples show that the use of any scalar part of Einstein's equations will not 
give a closed system of ordinary differential equations (see also Kofman & Pogosyan 1995). 
For the averaged variables we can also not expect this in view of Corollary 1: it states the 
equivalence of the averaged dynamics to the dynamics for the local field quantities. If this 
could be achieved in full generality, then this would be equivalent to having solved the full 
Einstein dynamics for the scalar parts, since the size of the domains could be arbitrarily 
chosen. 



3.3. Expansion law for closed universe models ? 

As discussed above an effort beyond the scope of ordinary differential equations for scalars 
is needed to close the system of equations for the average dynamics on any spatial domain. 

Notwithstanding, such an effort can be successful, if further constraints on the average 
dynamics are considered, most notably integral constraints which restrict the morpholog- 
ical characteristics of domains. Among them the integrated curvatures and, in particular, 
topological constraints that arise by restricting the Euler-characteristic of the surfaces 
bounding the domains. We already implied the topological constraint that the domains 
over which we average have to be simply-connected. In Newtonian cosmology (Buchert 
& Ehlers 1997) we have established a global criterion: if we extend this simply-connected 
domain to the whole Universe having topologically closed space sections (e.g., toroidal 
models), then this results in Qj) = on the closure scale. It is therefore to be expected 
that such a constraint may also close the system of averaged equations in general rela- 
tivity. We do not necessarily have the simple Newtonian condition. Corollary 3 states a 
general connection between Qx> and {TZ)xi'- a vanishing 'backreaction' would imply that all 
the contributions of the local curvatures that are produced by the inhomogeneities obey 
the "conspiracy" to sum up to the standard value oc a(t)~^ (where a{t) is a solution of a 
standard FRW cosmology). 

In curved spacetimes it is not straightforward to establish such a constraint and the 
line of arguments given in the Newtonian treatment is not conclusive in the present context. 
To illustrate this statement let us look at the extrinsic curvature tensor according to its 
definition as the (4— dimensional) covariant spatial derivative of the 4-velocity. Invariants 
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built from Kij and, consequently, the expression Qx>, cannot be written as total covariant 
divergences of vector fields in the hyper surfaces. As an example we look at the trace of 
Kij , K = — w".Q, = —9; the value of K on the hypersurfaces is covariantly defined, but the 
vector field -u* vanishes according to our spacetime foliation. A similar problem arises in 
the case of the second invariant. 

Forthcoming efforts should be directed towards finding a topological closure condi- 
tion for the hypersurfaces in order to determine the global average properties of the world 
models. This problem is more involved, since we cannot expect that the domains remain 
simply-connected. Working in a 4-dimensional tube of spacetime that is bounded by 
space-like hypersurfaces and considering the limit of vanishing distance between these hy- 
persurfaces, Yodzis (1974) attempted to derive average properties of closed spa,ce sections. 
His argument is reviewed in Appendix C, where it is shown that topological restrictions 
do not enter and his result holds for arbitrary compact and simply-connected domains. 

We conclude: 

3. We were not able to produce an argument analoguous to the Newtonian treatment 
stating that the 'backreaction term' vanishes for topologically closed space sections, if 
integrated over the whole space. Without such an argument averaged inhomogeneous 
cosmologies cannot be identified with the standard FRW cosmologies on any spatial scale. 
To justify this identification as an approximation there is presently no sufficiently general 
quantitative result as to whether the 'backreaction' term could be neglected on some scale 
or, in words suggested by Corollary 3, whether the averaged curvature decouples from the 
inhomogeneities . 



3.4. Remarks on perturbation theory 

I stated above that until present we don't know any quantitative calculation which may 
justify neglection of the 'backreaction' term on some scale. The reader may object that 
there exist several approximate calculations of the 'backreaction' effect in perturbation 
theory. However, there are severe obstacles for perturbative calculations which we are 
going to discuss now. 

As an example I would like to comment on a recent calculation by Russ et al. (1997; 
see also the references therein): based on the system of equations (10a), (10b) (using a 
rescaling of the metric and a split into background and perturbations) the 'backreaction' 
term was calculated within a second-order perturbation approach. In order to make the 
calculations concrete, Russ et al. have assumed periodic boundary conditions on some (very 
large) domain. Looking at their expression for the 'backreaction' term (B4) it is evident 
from their Eqs. (B7) and (B8) that 'backreaction' (in the sense defined in the present 
paper) has to vanish identically: together with mass conservation (their Eq. (BIO)) the 
introduction of periodic boundary conditions already leads to the result that the scale 
factor ai) obeys the standard Friedmann equations. (It should be noticed that Qj) = is 
already sufficient to have ax)(t) = a{t).) Hence, according to Corollary 3, it is no surprise 
that the average Ricci scalar has to obey the expansion law of models with spatial Ricci 
scalar oc a~^ in any consistent treatment of periodic perturbations on spatially fiat space 
sections. Therefore, Eq. (3.1) of Russ et al. (loc. cit.) cannot give any quantitative result 
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about the global value of the 'backreaction' term: it vanishes by assumption and the scale 
factor is given by the standard FRW cosmologies. Note also that the introduction of a 
Fourier transform, or a decomposition into plane waves, respectively, is only meaningful in 
the case of spatially flat space sections, i.e., also the averaged Ricci scalar has to vanish. 
It must be noted that the notion of 'backreaction' as implied by Russ et al. (loc. cit.) 
is slightly different from that in the present work. Any departure from the flat FRW 
cosmology in curvature (quantified by perturbation theory) may already be interpreted as 
'backreaction' (Kasai, priv. comm.). 

This attempt illustrates the possibly cyclic nature of calculations of the 'backreaction' 
term: if we start with spatially flat space sections and a model for the inhomogeneous 
deformation tensor, the standard methods of treating the perturbations as periodic on 
some scale already constrain the cosmology to one without 'backreaction' (measuring the 
deviations from the family of FRW cosmologies). Note that in perturbation theory the 
flrst-order perturbations are sources of higher-order perturbations and, as demonstrated 
by Russ et al. (loc. cit.), a large class of periodic flrst-order and, in turn, a large class of 
second-order perturbations on a flat hypertorus average Qj) to zero. 

In a realistic situation the domains on which one averages are not on average Ricci flat. 
Large domains (e.g. of the size of the Universe) may not be easily compactifled to make 
global statements about the evolution of the Universe: non-trivial topological spaceforms 
have to be considered and simple periodic boundary conditions are no longer useful. For 
some further remarks see (Buchert 1997). 



3.5. Global structure versus local models 

The insight gained from the set of generalized Friedmann equations (10) may be focussed 
in two ways: flrst, we are interested in the global structure of inhomogeneous cosmologies 
and, second, we would like to know more about average properties of spatial portions of 
the Universe without severely restricting the dynamical model. 

As for the flrst point, a globally non-vanishing 'backreaction' that may be small at early 
epochs of the Universe' evolution could, on the scale of the size of the Universe, slowly build 
up due to an ampliflcation of inhomogeneities. The more large-scale structure develops, 
the more the whole Universe might undergo global changes in morphology including the 
possibility of topology change. In order to analyze global changes during inflationary 
stages the present matter model 'dust' has to be generalized (which is the subject of a 
forthcoming work. Paper II). At later epochs constraints from the microwave background 
anisotropies can be used to give upper limits on the 'backreaction' characteristic Oq (Eq. 
12d): if we accept that the microwave background dipol is only due to our proper motion 
against an isotropic CMB reference frame, then limits on the global shear parameter 



icr') 



CMB /.p N, 

- — ^rr2 (15a) 

may be related by the assumption that, on the CMB scale, 

~ (^)cMB ; (156) 
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then, on the this scale, 

Oq ^ E^ . (15c) 

Maartens et al. (1995) have given upper hmits on the shear parameter for a Bianchi-type 
universe, in which case E < 10~^ is a tight constraint on the global magnitude of the 
'backreaction' characteristic at the epoch of last scattering (see also: Wainwright & Ellis 
1997). As a disclaimer we note that an average over inhomogeneous models is performed 
in Eq. (15a), and the average model must not necessarily isotropize as the Bianchi-type 
models (except type IX) do. Therefore, care must be taken in using such constraints at 
times after last scattering. 

As for the second point, the expansion laws discussed in Appendix B and similar 
relations (calculated, e.g., from perturbation theory, or hybrid models employing pertur- 
bation theory on large scales, but including the full nonlinearities on small scales, Takada 
& Futamase 1999) provide a more general architecture for the study of hierarchical cos- 
mologies, understood in the sense of models which do not single out the non-generic case 
of scale-independent mean density, as the standard model does. The focus here is on the 
effective dynamics of portions of the Universe on some spatial scale including the possi- 
bility of statistically studying ensembles of spatial domains (see: Buchert et al. 1999 for 
an investigation within Newtonian cosmology). With regard to the old ideas of hierarchi- 
cal cosmologies the suggestion by Wertz (1971) may be put into perspective. It relies on 
spherically symmetric domains which depend on their own parameters of a standard FRW 
cosmology. In this line the expansion laws (B.3) (and more general relations) can also be 
applied to finite domains having their individual parameters. These parameters belong to 
generalized FRW cosmologies that include averages over inhomogeneities encoded in an 
additional 'backreaction' parameter (Eq. (12d)). The relevance of the characteristics (12) 
on a finite domain for the interpretation of volume-limited observational data, where we 
cannot a priori assume that the surveyed volume is a portion of a standard Hubble fiow, 
is obvious. 

A similar view applies to so-called collapse models like the spherical "top-hat" model 
(e.g., Peebles 1980): on smaller spatial scales expansion fluctuations may become dominant 
in an overdense domain leading to collapse. The transition when Qx> moves through zero as 
we come from large scales can be used to mark the scale of "decoupling of inhomogeneities 
from the global expansion" . In this context 'backreaction' models provide straightforward 
generalizations of the top-hat model. While Birkhoff's theorem lies at the basis of the 
spherical model, the average dynamics including 'backreaction' is not restricted by sym- 
metry assumptions. Thus, the averaged equations furnish a general framework with which 
one can describe the effective dynamics of individual collapsing or expanding domains. 
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Appendix A: NeAvtonian Analogues 



The general expansion law in Newtonian cosmology reads (Buchert & Ehlers 1997): 

3^ + 4.G^-A=S. . (-4.1) 

As in the main text, Md denotes the total (conserved) rest mass contained in V, and Qj) 
is the same expression as Eq. (13c). This equation is identical to (10a), as we already 
pointed out in (Buchert, Ehlers 1997). 

The integral of Equation (A.l) with respect to time can be performed and yields the 
generalized form of Friedmann's differential equation for the first derivative of the scale 
factor (Buchert 1996 - with a different sign convention for Qx))'- 

a|, + k-v SttGMd A _ 1 /"* , / ^ '^2 



ttv ^V-D^a^ 3 3a|, 



J dt' Qv ^o^lit') , (A.2) 



where kx> is a (domain-dependent) integration constant. 

Comparing with the general relativistic equation (10b) we find the analogy: 

|-3i/**'2-i!74(*')4(W. + S.) . (-4.3) 

The time-derivative of Eq. (A3) is equivalent to the integrability condition (13) in Corol- 
lary 3. We may view Eq. (A3) as an integral of (13b). Eliminating the average curvature 
from this integral, Eq. (A3), and inserting it into the integrability condition (13b) formally 
results in a differential equation for Qx> alone, which, however, reduces to an identity. 

Notice that we cannot separately identify the integration constant kx> with the average 
Ricci scalar, since this would determine the evolution of the average curvature and the 
'backreaction' term to the special solution (B.3). We might be able to show that, e.g., 
the solution (B.3) could also be found within the Newtonian framework for a special type 
of deformation of the domain's boundary. However, we cannot conclude that for the 
subclass of Newtonian cosmologies, which can be obtained from the Newtonian limit of 
the corresponding class of GR solutions, the solutions (B.3) would be the correct limit; the 
limit (c ^ 00; 7^ ^ 0) is not obvious in the expression c^lZ. 
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Appendix B: Examples of Expansion Law^s 

From Equation (13b) we conclude that, 

first, for on average spatially fiat domains the 'backreaction' can be integrated to give 

Q^(t) = Q0,a^6 ; Q'^:=Qv{to) • (B.la) 

Using the integral (B.la) we write down a closed equation for the scale factor in this case: 






3^+47rG-^-A = ^ . [B.lh) 



In view of Eq. (10b) the integral of (B.lb) is given by 






3 (nil) _87rG'-^-A = -^ . (S.lc) 



Hence, the problem is reduced to a quadratur. 

Second, for (on some spatial domain) vanishing 'backreaction' we obtain the special case 
of conformally constant curvature models; the average scalar curvature is inversely pro- 
portional to the square of the "radius of curvature" on the domain, 

(7^)x)(t) = (7^)^ a^2 ; (7^)0, := (7^)x)(to) , {B.2a) 

where aTi{t) = a{t) is a solution of the standard FRW models. We are faced with the situ- 
ation that the domain represents on average a small FRW universe with its own domain- 
dependent parameters. However, here, this result holds for any spatial domain on which 
the 'backreaction' vanishes, i.e. inhomogeneities are present and their fiuctuations can 
even be large, but they have to compensate each other. 

Third, one obvious solution of Eq. (13b) in the case of non-vanishing average scalar 
curvature and non-vanishing 'backreaction' is given by 

Qvit) = Qla^'' ; (7^)p(t) = (7^)°, a^2 . (S.3a,6) 

here, (7^)i? is proportional to an "effective radius of curvature" of the domain, and 
dvit) 7^ a{t). This solution features independent evolution of average Ricci scalar and 
'backreaction': the spatial domain behaves like a small "almost FRW" universe, still being 
characterized by its own parameters, which are exclusively determined by the values of the 
fields inside the domain. Although special, this case offers the possibility of understanding 
some properties of the 'backreaction' effect. 

At first glance, it might look wrong that the dynamics of any patch of matter can be 
described independently of the environment; the non-local gravitational infiuence from the 
matter outside the domain seems not to have impact on the dynamics of the domain. This 
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interpretation is, however, misleading: although we have to specify only initial data within 
the domain in, e.g., the solutions to (B.3), we still have to solve the constraint equations 
for these initial data which is a non-local operation and involves also the fields outside 
the domain under consideration. Still, solutions to (B.3) uniquely describe the averaged 
dynamics of the domain for all times, once the initial data are specified, without solving 
the constraints at later times. 

All of our examples restrict the generality: solution (B.l) is found by the require- 
ment of on the domain vanishing average Ricci scalar. Hence, in view of the Hamiltonian 
constraint (4a) the inhomogeneities have to obey a special relation between the rest mass 
density, the cosmological constant and the second scalar invariant of the expansion tensor. 
Eqs. (B.3) together with the generalized Friedmann equations determine a more general 
class of motions. Still, in general, we expect that the evolution of 'backreaction' is coupled 
to the evolution of the averaged spatial curvature in a more complex way. Indeed, as Eq. 
(13a) shows, even for initially vanishing average curvature, there is generation of curvature 
in the course of structure formation, since amplification of inhomogeneities builds up the 
'backreaction' term. 

We finish the discussion of special expansion laws by giving a useful formula for the 
dynamical relation of the average characteristics (12). Combining (12) with Corollary 3 
and using the Theorem we obtain: 

Oq + &HTyQQ{l - Ofc - Qq) + Ofc + 2if:D0fc(l - Ofc - Qq) 

-3Hr){l-nA-^k-^Q){^k + ^Q) = . (BA) 

One example may express a warning that the average characteristics in the present case 
are dynamically tightly related and should not be treated independently: let us assume 
that the cosmological term vanishes (which remains an independent parameter) , and that 
the restmass density characteristic remains constant in time mimicking the situation in 
an Einstein-de Sitter universe. Then, in the simplest case flrn = 1? the curvature and 
'backreaction' characteristics have to compensate each other, Ofc + Oq = 0, and Eq. (B.4) 
reduces to 

Oq + QHt^VLq + Ofc + 2H'Dnk = . 

Eliminating one of the characteristics from this equation shows that the other has to vanish 
identically, reducing the average model to the standard Einstein-de Sitter universe. There 
exists no generic inhomogeneous cosmology with O^ = 1 throughout its evolution. 
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Appendix C: Yodzis' Argument 

Let us view the averaging problem within a 4-diniensional tube of spacetime that is swept 
out by a compact 3-dimensional domain V and bounded from above and below by space- 
like hypersurfaces ti = const., t2 = const. From the Ricci identity 

we obtain after contraction in a and /3 and multiplication with u"^: 



u u .p — w u 



P,,"- A — _P .,,«^,/^ -l_ ^,7 ^2_^,7 ,,S 



^];« = -i?«/3t^"«^ + «^)'-«V;7 ■ (^-1^) 



Recalling the definition of the second invariant of Kij (8b), we rewrite Eq. (C.lb): 

211 = K-u^.^ - uf^W'.p].^ + Rc,(3U°'u'^ . (C.lc) 

After averaging we obtain for the 'backreaction' term: 



2 



J JcfX [u'^uf'.p - ufu'^.p],^ ~W^{j "^"^'^ ""^'/s) + {Ro.i3u''u'')v . (C.2) 



We now need to evaluate the total 4-divergences in Eq. (C.2) on the hypersurfaces. Yodzis 
(1974) gave the answer using the following argument: perform the volume integral over the 
4-dimensional tube of spacetime F and let then the distance between the hypersurfaces 
tend to zero, ^2 — ti = e ^ 0. Applying Green's theorem on the integrals of the invariants. 



a 
UL U 

r(ti,t2) 



dV u',^ 



- f Jd^X u'^uc, + / Jd^X u'^uc, = Vvit2) - Vviti) ; (C.Sa) 

[ dT[u"u^.0-uf^u"s]-c = 

- f Jd^X [u'^u^.p - -u^-u^.^jwa + / Jd^X [u'^u'^.fj - uf^u''.p]uc, , (C.36) 

he arrived at the following result by deviding Eqs. (C.3a,b) by e and taking the limit 
e ^ (in our notations and conventions): first, he obtains the familiar equation (9d) for 
the average of the first invariant 1 = 9: 

Jd^X e = {e)vVv ; (C.3c) 

V 
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for the average of the second invariant II he derives: 

/ Jd^X [-R^pu'^u'^ + 211] = {{9)t)Vv)' . {CM) 

Jv 

Although Yodzis (loc. cit.) seems to imply that this result only holds for closed 3-spaces, we 
can immediately see that these equations and especially Eq. (C.3d) hold for any compact 
and simply-connected domain in the hypersurfaces: from Einstein's equations (la) we 
have Rapu^u^ = AnGg — A. Thus, Eq. (C.3d) is equivalent to the averaged Raychaudhuri 
equation in Corollary i, Eq. (lie). Calculating Qx> from Eq. (C.2) we accordingly obtain 
Eq. (10a) of the main text. 
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